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ABSTRACT 

Cyanogen (NCCN) is the simplest member of the series of dicyanopolyynes. It has been hypothesized that this family of molecules 
can be important constituents of interstellar and circumstellar media, although the lack of a permanent electric dipole moment prevents 
its detection through radioastronomical techniques. Here we present the first solid evidence of the presence of cyanogen in interstellar 
clouds through the detection of its protonated form toward the cold dark clouds TMC-1 and L483. Protonated cyanogen (NCCNH + ) 
has been identified through the 7 = 5 — 4 and 7 = 10-9 rotational transitions using the 40m radiotelescope of Yebes and the IRAM 
30m telescope. We derive beam averaged column densities for NCCNH 4 of ( 8.6 ± 4.4) x 10 10 cnT 2 in TMC-1 and (3.9 ± 1.8) x 
10 10 cm ” 2 in L483, which translate to fairly low fractional abundances relative to H 2 , in the range (1-10) x 10” 12 . The chemistry of 
protonated molecules in dark clouds is discussed, and it is found that, in general terms, the abundance ratio between the protonated 
and non protonated forms of a molecule increases with increasing proton affinity. Our chemical model predicts an abundance ratio 
NCCNH + /NCCN of ~ 10” 4 , which implies that the abundance of cyanogen in dark clouds could be as high as (1-10) x 10 ” 8 relative 
to H 2 , i.e., comparable to that of other abundant nitriles such as HCN, HNC, and HC 3 N. 
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1. Introduction 


Nitriles, i.e., molecules containing a functional group —C=N, 
are present in diverse astronomical environments. In particular, 
cyanopolyynes, H-(C=C)„-C=N, are commonly found in cold 
interstellar clouds (IBroten et al.H 978: IBell et akl l997). circum¬ 
stellar envelopes around evolved stars, especially in carbon-rich 
objects (Winnewisser & Walmsley 1978: lPardo et al. 1 2005). and 
planetary atmospheres with a high content of carbon and nitro¬ 
gen, as that of Titan (IKunde e t al. 1981: ICoustenis et al.1 1991). 


It has been suggested that dicyanopolyynes, molecules con¬ 
taining two cyano groups, N=C-(C=C)„—C=N, could be abun¬ 
dant in interstellar and circumstellar clouds (Kolos & Gra bowskil 
2000; iPetrie et akl 2003). However, these molecules cannot be 
detected through their rotational spectrum because they lack 
a permanent electric dipole moment. The simplest member of 
this series, cyanogen (NCCN), is thought to be a major pre- 
the CN radical observed in cometary comae (e.g., 
2005), although its detection in a comet still remains 
challenging. Moreover, infrared observations carried out with 
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* Based on observations carried out with the IRAM 30m Telescope 
and the Yebes 40m Telescope. IRAM is supported by INSU/CNRS 
(France), MPG (Germany) and IGN (Spain). The 40m radiotelescope 
at Yebes Observatory is operated by the Spanish National Geographic 
Institute (IGN. Ministerio de Fomento). 


the Voyager 1 have identified NCCN in the atmosphere of Titan 
(IKunde et al.l 1981) and the larger homologue NC 4 N has long 
been thought to be present as well (Jolly et al. 2015). It is also 
worth noting that a chemical cousin of cyanogen in which one N 
atom is substituted by a P atom, NCCP, has been tentatively iden¬ 
tified in the C-rich envelope IRC +10216 ( [Agundez et al.| 2014). 

Given the difficulty to directly detect NCCN and larger di¬ 
cyanopolyynes in cold interstellar clouds, it has been proposed 
that indirect strategies to probe their presence would be to ob¬ 
serve chemically related molecules such as the polar metastable 
isomer CNCN or the protonated form NCCNH + ( Petrie et~ai 7 
2003). Here we report the first detection in the interstellar 
medium of NCCNH + , the protonated form of cyanogen, toward 
the cold dark clouds TMC-1 and L483. 


2. Observations 

Protonated cyanogen is a highly polar linear molecular cation. 
Its rotational spectrum has been characterized in the laboratory 
from microwave to millimeter wavelengths 
1991: [Gottlieb et akl 2000) and its electric 
been calculated as 6.448 Debye (IBotschwina & Sebaldl 1990). 
We have observed the 7 = 5-4 transition at 44.4 GHz with the 
Yebes 40m telescope and the 7 = 10-9 transition at 88.8 GHz 
with the IRAM 30m telescope toward the cold dark clouds 
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Table 1 . Observed line parameters of NCCNH + in TMC-1 and 
L483 


Transition 

Frequency 

(MHz) 

Vlsr 
(km s~‘) 

Av 

(km s~‘) 

J T l dv 

(K km s~') 

TMC-1 

7 = 5-4 

7 = 10-9 

44379.850 

88758.108 

+5.68(8) 

+5.91(9) 

0.72(13) 

0.48(12) 

0 .011(2) 

0.006(1) 

L483 

7 = 5-4 

7 = 10-9 

44379.850 

88758.108 

+5.11(5) 

+5.31(9) 

0.26(4) 

0.49(13) 

0.006(1) 

0.017(4) 


Numbers in parentheses are lcr uncertainties in units of the last digits. 


TMC-1 and L483 at the positions observed by Mar celino et al.l 
(2005, 2007) and |Agundez et al. (2008), respectively. 


2.1. I RAM 30m 

The IRAM 30m observations were carried out using the EMIR 
3 mm receiver and the fast Fourier Transform spectrometer with 
a spectral resolution of 50 kHz. The frequency switching tech¬ 
nique was used to optimize the telescope time. The half power 
beam width (HPBW) of the IRAM 30m telescope at 88.8 GHz 
is 2773. The observations of TMC-1 are part of a spectral line 
survey at 3 mm (IMarcelino et al.l 2005. 2007, 2009). Most of 
the observations at 88.8 GHz used in this article were taken 
in February 2012, when the good weather conditions resulted 
in a system temperature of ~ 70 K. More details are given in 
ICernicharo et al.l (2012). The observations of L483 were taken 
from September to November 2014 during an observational 
campaign aimed at observing negative ions in dense molecular 
clouds and are described in detail in Agundez et al. (2015). 


2.2. Yebes 40m 

The 40m radiotelescope is located at Yebes (Guadalajara, Spain) 
at 980 m above sea level. Weather conditions are typically dry; 
the ammount of precipitable water vapor ranges from 4 mm in 
winter to 14 mm in summer. The antenna has an homological 
design and its optics is that of a Nasmyth radiotelescope with 
several receivers in the Nasmyth cabin. The antenna is equipped 
with a cryogenic single pixel dual polarization 45 GHz receiver 
built at Yebes, whose instantaneous bandwidth ranges between 
41 and 49 GHz. The IF only processes bands of 500 MHz and/or 
1500 MHz and its signal is injected into a fast Fourier Transform 
spectrometer with several modules, which can be configured 
to provide 16384 channels on a bandwidth of 100 MHz or 
500 MHz. Observations reported in this article used two mod¬ 
ules 100 MHz wide with a spectral resolution of 6 kHz. Spectra 
were later on smoothed to a spectral resolution of 12 kHz. 

The observations toward TMC-1 and L483 were carried 
out in 9 and 12 periods, respectively, of 5-7 h during March 
and April 2015. System temperatures ranged between 120 and 
200 K. On-off observations were done with an integration time 
of 60 s and the off position located 600" away in right ascension. 
Calibration was repeated every 20 minutes and consisted of ob¬ 
serving the sky and a hot load at ambient temperature. The sky 
opacity was estimated using the weather conditions measured 
400 m away from the antenna and the atmospheric transmission 
model ATM (ICernicharol 1985; IPardo et al.l 2001). We believe 
that the sky opacity has a maximum error of 10 % (estimated by 
comparing with skydip measurements). To keep a good point¬ 
ing and focus, continuum observations were performed every 20 




Fig. 1 . Spectra of TMC-1 and L483 showing the emission lines 
assigned to the J - 5-4 and J = 10-9 rotational transitions 
of NCCNH+, lying at 44.4 GHz and 88.8 GHz, respectively. The 
T* a rms of the spectra at 44.4 GHz is 0.0045 K per 12 kHz chan¬ 
nel for both sources, while at 88.8 GHz the T* A rms is 0.0029 K 
and 0.0088 K per 50 kHz channel, for TMC-1 and L483, re¬ 
spectively. The antenna temperature scale can be converted to 
main beam brightness temperature by dividing by (B clt / F ci (), 
which takes a value of 0.36/0.90 for the Yebes 40m telescope 
at 44.4 GHz and of 0.82/0.90 for the IRAM 30m telescope at 
88.8 GHz. 


minutes towards UOri (for TMC-1) and OH26.5+0.6 (for L483). 
Pointing and focus are estimated to be accurate within 5". The 
HPBW of the Yebes 40m telescope at 44.4 GHz is 4276. 

The final spectra of TMC-1 and L483 were obtained by av¬ 
eraging 2440 and 2338 individual spectra, respectively, half of 
them with left(right) circular polarization. The intensity scale is 
given in antenna temperature (T*), which corrects for the sky 
opacity and forward efficiency (90 %). The main beam efficiency 
has been estimated to be 36 % in an elevation range between 20° 
and 80° from observations toward Venus and Saturn. Calibration 
can be considered correct within a 10 % of uncertainty, mostly 
coming from the aperture efficiency and atmospheric conditions. 


3. Results 

The emission lines observed toward TMC-1 and L483 and as¬ 
signed to the 7 = 5-4 and J = 10-9 transitions of protonated 
cyanogen are shown in Fig. |T] Line parameters derived from 
Gaussian fits using GILDAS are listed in Table |T| The obser¬ 
vations indicate that in TMC-1 the 7 = 5-4 line is more intense 
than the 7 = 10-9, while in L483 the contrary is found. This fact 
suggests that NCCNH + has probably a higher rotational temper¬ 
ature and/or a more compact spatial distribution in L483 than in 
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TMC-1. The / = 5 -4 and J = 10-9 transitions have upper level 
energies of 6.4 K and 23.4 K, respectively, and are observed with 
main beam sizes of 4276 and 27" 3, respectively. 

Since we do not have information on the spatial distribution 
of NCCNH + in these two sources, we may assume that the emis¬ 
sion size fills the main beam of the IRAM 30m and Yebes 40m 
telescopes to derive a first order estimate of the column densities 
of protonated cyanogen. In TMC-1, we derive a beam-averaged 
column density of (8.6±4.4)x 10 1H cm -2 and a rotational temper¬ 
ature of 5.8 ± 0.8 K. In L483, we derive a somewhat smaller col¬ 
umn density, (3.9+1.8) x 10 1() cm' 2 , and a higher rotational tem¬ 
perature, 13.4+4.5 K, which however has an important degree of 
uncertainty. If we adopt H 2 column densities of 1 x 10 22 cm' 2 in 
TMC-1 (tCernicharo & Guelinl 19871 and 3 x 10 22 cm' 2 in L483 
(ITafalla et al.l 2000; see also [Agundez et alT| 2015), we end up 
with fairly low fractional abundances relative to IT, 8.6 x 10 12 
in TMC-1 and 1.3 x 10 -12 in L483. 

4. Discussion 

The detection of NCCNH + is a good indication of the presence 
of NCCN in dark clouds. However, to have an idea of how abun¬ 
dant is cyanogen in this type of sources it is necessary to have 
a look at the chemistry of protonated molecules in general, and 
NCCNH + in particular. In a simplified chemical scheme, a pro¬ 
tonated molecule MH + can be formed by proton transfer to M 


XH + + M —> MH + + X, 


(1) 


MH + + e —> products, 


( 2 ) 


where the products can be assorted neutral fragments. Within 
this simpe chemical scheme, at steady state we have 


[MH + ] k PT [XH + ] 


[M] 


K DR 


[e ] 


(3) 


where k PT and k DR are the rate constants of the reactions of pro¬ 
ton transfer and dissociative recombination, respectively. Eq. 0 
suggests that high proton affinities of M would tend to enhance 
the rate of proton transfer (by increasing the number of avail¬ 
able proton donors XH + and probably the associated rate con¬ 
stants kpj) and so the abundance ratio [MH + ]/[M], The above 
chemical scheme may however be not suitable for all protonated 
molecules, making necessary to build a chemical model that in¬ 
cludes all relevant reactions to provide more precise predictions. 

We have built a pseudo time-dependent gas phase chemical 
model of a dark cloud adopting standard physical parameters 
(T k = 10K, n H = 2x 10 4 cm” 3 , £ = 1.3xl0 -17 s~ l ,A v = 30) and 


“low metal” elemental abundances (see Agundez & Wakelam 
2013). We have adopted the UMIST RATE12 reaction net- 


calculated c = 5 klO 17 s 1 
calculated £ = 1.3 ?<10 17 s 1 


cA 


HCN/HNC 


£ £ 


i nh. 


NCCN 

£ 


where M must have a higher proton affinity than X in order to 
make the reaction exothermic. On the other hand, in cold dark 
clouds molecular cations are usually depleted by dissociative re¬ 
combination with electrons 


work (Vic Elro y et al.| 2013) w ith a subset of reactions involv¬ 
ing HCCN from iLoison et al.l (2015). We find that the chemi¬ 
cal model validates the leading role of dissociative recombina¬ 
tion, reaction 0, as the major destruction process of protonated 
molecules, except for N 2 H + and HCOj, which are also depleted 
to an important extent by proton transfer to CO. We also find 
that, depending on whether the proton affinity of M is below or 
above that of CO, the main proton donor XH + in reaction 0 
is either H i or HCO + , respectively. However, some protonated 
molecules such as HCO,/, HCS + , and HCAS 1 are mainly formed 


120 140 160 180 200 

proton affinity of M (kcal mol -1 ) 

Fig. 2. Abundance ratios [MH + ]/[M] between the protonated and 
non protonated form of some molecules as a function of the pro¬ 
ton affinity of M. Calculated values at steady state (reached after 
some 10 s yr; see Fig.[3j), adopting two different cosmic-ray ion¬ 
ization rates are compared with values derived from observations 
of TMC-1 ( [Agundez & Wakelam 2013). In the case of NHJ, the 
observed value is derived from observations of the monodeuter- 
ated species in Barnard 1 (ICernicharo et al.l 20131. The values 
for HCN/HNC include both isomers and are plotted at the mean 
value of the proton affinities of HCN and HNC, 170.4 kcal mol 1 
and 184.6 kcal mol -1 , respectively. Note the trend of increasing 
[MH + ]/[M] with increasing proton affinity of M. 


by ion-neutral processes other than reaction (jTJ, in which case 
Eq. 0 underestimates the abundance ratio [MH + ]/[M]. 

Appart from NCCNH + and the widespread ions HCO + 
and N 2 H', some other few protonated molecules have been 
observed in cold dark clouds: HCS + ( Thaddeus et al . 1981), 
HCNH+ (ISchilke et al.l 1991), HC 3 NH + (|Kawaguchi et al. 
1994), HCO+ ([Turneret al.l l999: ISakai etaP 20081. andNH,D + 
(ICernicharo et al.1 20131. The abundance ratio [MH + ]/[M], when 
available from observations, is well reproduced by the chemical 
model for M = CO and NH 3 , although it is underestimated for M 
= HCN/HNC, HC 3 N, and especially CS (see Fig. [2]). Taking into 
account that HCNH + , HC 3 NH + , and HCS + are mainly destroyed 
by dissociative recombination with electrons, whose rate con¬ 
stants are well known from experiments (ISemaniak et al.1 2001: 
Geppert et al. 2004; Montaigne et al. 2005), and that the rate 
constants of formation by proton transfer from HCO + are usually 
well constrained experimentally (IAnicichl 2003). it is likely that 
the underestimation occurs because the chemical model misses 
important formation routes to HCNH + , HC 3 NH + , and HCS + . 

The abundance ratios [MH + ]/[M] are sensitive to the degree 
of ionization, and thus, to various physical parameters such as 
the cosmic-ray ionization rate £ (the sensitivity is more marked 
for low proton affinities; see Fig. [2]) and the volume density of 
particles (the higher the density, the lower the ionization fraction 
and thus the lower the importance of protonated molecules). It is 
interesting to note that both the chemical model and the obser¬ 
vations suggest a trend in which the abundance ratio [MH + ]/[M] 
increases with increasing proton affinity of M. Since the destruc¬ 
tion of MH + is controlled by dissociative recombination with 
electrons, whose rate constants are usually similar within one 
order of magnitude, the trend can be explained in terms of an 
enhanced formation rate of MH + with increasing proton affin¬ 
ity of M, resulting from the appearance of multiple formation 
pathways to MH + through exothermic ion-neutral reactions of 
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time (yr) 

Fig. 3. Calculated abundances of NCCNH + and NCCN as a 
function of time using the UMIST RATE 12 chemical network 
(solid lines) and after removing the reactions of destruction of 
NCCN with C and N atoms (dashed lines), in both cases adopt¬ 
ing £ = 1.3 X ICC 17 s -1 . The range of abundances of NCCNH + 
observed in TMC-1 and L483 is also indicated. The calcu¬ 
lated abundance ratio [NCCNH + ]/[NCCN] is shown in the upper 
panel as a function of time. 


seems the most immediate step to better constrain the chemistry 
of cyanogen in cold interstellar clouds. Additional uncertain¬ 
ties could come from other reactions that affect the precursors 
HNC and CN. For example, the systematic use of enhanced rate 
constants for most ion-polar reactions, adopted in the UMIST 
RATE 12 network, results in lower abundances for HNC, NCCN, 
and NCCNH + , as compared with the use of a chemical network 
in which these rate constant enhancements are not adopted (see, 
e.g JWoodall et al.l 2007). Formation of cyanogen on the surface 
of dust grains followed by some non-thermal desorption process 
is also a possibility to explore. 

If we trust the abundance ratio [NCCNH + ]/[NCCN] of ~ 
10 4 calculated by the chemical model, the abundance of NCCN 
in dark clouds could be as high as (1-10) x 1CL 8 relative to H 3 , 
i.e., comparable to that of other abundant nitriles such as HCN, 
HNC, and HC 3 N. 
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proton transfer or other type. We may therefore expect to find 
other abundant protonated molecules MH + in dark clouds, as 
long as M has a high proton affinity and a sufficiently high 
abundance. Molecules such as H 2 COH + (previously detected 
by lOhishi et al.i (1996) toward warm, but not cold, clouds), 
OCSH + , and CnHj are potentially detectable. However, the most 
promising candidate is probably HC 2 S + , given that it has a high 
proton affinity and dipole moment (Barrientos & Largo 1991; 
Maclagan & Sudkeaw 1992: lPuzzarinil 2008) and it is predicted 
to be just a few times less abundant than C 2 S, whose column 
density in dark clouds is in the range (1-100) x 10 12 cirT 2 
(IFuente et al.l l990: ISuzuki et al.1 1992). 

The chemical model predicts a low abundance of NCCNH + 
(see Fig. [3]), although it is not clear whether it is because the 
model underestimates the abundance of NCCN or the abundance 
ratio [NCCNH + ]/[NCCN]. According to the chemical model, 
NCCNH + is formed by proton transfer from HCO + to NCCN 
and destroyed by dissociative recombination with electrons, the 
rate constants of which are not known, although they are un¬ 
likely to be radically different from the values guessed in the 
UMIST RATE 12 reaction network. We thus expect the abun¬ 
dance ratio [NCCNH + ]/[NCCN] to be around 10 4 , as predicted 
by the chemical model, or higher, if important formation routes 
to NCCNH + are missing in the model. If the abundance ratio 
[NCCNH + ]/[NCCN] is correctly predicted, then the low frac¬ 
tional abundance predicted for NCCNH + must arise from an 
underestimation of the abundance of NCCN. In the chemical 
model, cyanogen is essentially formed by the reaction 


HNC + CN —> NCCN + H, 


(4) 


for which iPetrie et aht (2003) estimate a rate constant of 2 x 
s -1 , and depleted by reaction with C and N atoms. 
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whose rate constants at room temperature are inferred to be 


a few 10 11 cm 3 s 1 (Whyte & Phillips 1983; |Safrany & J aster 


1968). The calculated abundance of NCCN (and NCCNH+) is 
very sensitive to these rate constants (see Fig. [3ji. Therefore, a 
better understanding of the low temperature chemical kinetics 
of reaction 0 and the reactions of NCCN with C and N atoms 
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